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Abstract: In this paper we solve the equation for nucleus- nucleus scattering in the BFKL Pomeron 
calculus, suggested by Braun [1]. We find these solutions analytically at high energies as well as numerically 
in the entire region of energies inside the saturation region. The semi-classical approximation is used to 
^ pelect out the infinite set of the parasite solutions. The nucleus- nucleus cross sections at high energy are 
• 'H Estimated and compared with the Glauber-Gribov approach. It turns out that the exact formula gives 
''jjj the estimates that are very close to the ones based on Glauber-Gribov formula which is important for the 
_ practical applications. 
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1. Introduction 

The goal of this paper is to find the solution to the equations for nucleus-nucleus collision that have been 
derived in Ref. [1]. We continue the attempts, taken in Refs. [2-5], to study these equations and to search 
the general method of solving them. 

Nucleus-nucleus scattering gives the most informative example of dense - dense parton system inter- 
actions in which we can see the main prediction of Color Glass Condensate/saturation approach [7-12]. 
However, in spite of the fact that we know quite well qualitative features of nucleus-nucleus scattering 
(see Refs. [13]) CGC/saturation approach suffers by the absence of the evolution equation that gives us 
a possibility to find the scattering amplitude at high energy. On the other hand we know quite well the 
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initial condition for such an evolution [9]. Fortunately, the second approach to the high energy QCD; the 
BFKL [14, 15] Pomeron calculus, gives the equations for nucleus-nucleus scattering at high energy. For 
dilutc-dcnsc parton system scattering both approaches: BFKL Pomeron calculus and Color Glass Conden- 
sate (CGC), lead to the same nonlinear equations [1,16]. Therefore, we can hope that the equations given 
by BFKL Pomeron calculus would be proven in the framework of CGC. 

In the next section we give the brief review of the equations derived in Ref. [1] and discuss their 
main properties. This section does not contain any new results except section 2.3, and it is written for 
the completeness of presentation. In section 2.3 we consider the asymptotic solution to the problem at 
large values of rapidity Y in the framework of the semiclassical approach that has been developed by us in 
Refs. [4,5]. In the next section we show that the number of possible solutions has to be reduce to unique 
solution which is discussed in section 5. In section 6 we derive that nucleus-nucleus amplitude for the 
solution given in section 5. In conclusions we summarize our results and compare our solution with the 
numerical solutions of Refs. [2,3]. Unfortunately, the main equations were proposed a decade ago but we 
have only had five papers devoted to a search of the solutions (see Refs. [2-5] and this paper). 

2. The BFKL Pomeron calculus for nucleus-nucleus interaction at high energy 
2.1 Equations for nucleus-nucleus scattering 

The most economic and elegant form the BFKL Pomeron calculus has in terms of the functional integral [1] 



where Sq describes free Pomerons, Sj corresponds to their mutual interaction while Se relates to the 
interaction with the external sources (target and projectile). 

We will write these actions in the momentum representation [2, 5] which is defined as 




with S = Sq + Si + Se 



(2.1) 



$(xi,x2,r') 



$t {xi,X2,Y') 




(2.2) 



(2.3) 



Sq takes the form 



So = 64 (27r) 




(2.4) 



where I = Ink^ and 



'H^{k,b,Y) = as< 



cPk' 



ecPk' 



^^{k,b,Y) \ 



(2.5) 
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where as = (Nc/tt) as {Nc is the number of colours and the running QCD coupling as = 1/ yPo In yk"^ /Aqfj^ 
and /3o = (33 — 2n/) /127r with n/ is the number of the fermions). 
The interaction term 5"/ can be written as follows [5]: 

Si = 16^^^ / dY'l AdH I cPk | cl>t [-k^Y') $t {-k^Y') (| + l)' {k^Y') 

+ $(_fc,6,y')$(_fc,6,y')(^| + iy ^$t(A,,6,y')| (2.6) 



The equations for nucleus-nucleus scattering have been derived from the averaging of the equations of 
motion for the action of Eq. 



5S \ ^ / 6S 



.S^{k,b,Y')/ \6^^k,b,Y') 
where averaging is understood as 



(2.7) 



Deriving the equation of motion we assume that 



^^k,Y';b)) = (U{k,Y';b))Y (2.9) 



2 \ / / , \\2 



$ (/t, y ; 6) $t (/t, y ; 6) ) = (A;,y';&) ) X (A;,y';6) 



These identities are proven in the case of nucleus-nucleus scattering within accuracy of about (see 
Refs. [1,4,11]). We need to find the relation between fields ^ {k,b,Y') and <I>^ {k,b,Y') and the scattering 
amplitude owing to the single BFKL Pomeron exchange for Sj = which we denote N {k,kQ;b,Y') (where 
k and ko is the final and initial transverse momenta at rapidity Y and Yq, respectively). 

Taking into account Eq. ( |2.9D one can see that the variation with respect to <I>^ {k, b, Y') leads to the 
following equation of motion 



6 {So + Si)/6<^Hk,b,Y') 



64(27r) 



d 



+ 16 



'27ra1 



d 



dl ^ j df KdY 



{2n)' { 2c^t b,Y')(- + l) j^^{k,b,Y') + (^ + l 



Hj^ {k,b, Y') 

2 Q2 



(2.10) 



5/2 



d 



dl 



-k,b,Y') 
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Using 5*0 we easily see that 

N{k,ko,b;Y;Yo) = 27t^ as {HI, b, Y -Yq)) (2.11) 



For understanding the relation between field <I>^ and N {k, KQ,b;Y;Y()) where Kq is the transverse 
momentum at rapidity Y we use the equation [7, 17] 



NiKo,ko,b;Y;Yo) = J d%'dl (^Ko,k,b - b' ,Y - Y'^ N (k,ko,b' ,Y' - Yq^ (2.12) 

Eq. (|]l|) has more general meaning than for exchange of one Pomeron (see Ref. [4] for proof in the case 
of nucleus-nucleus scattering) : it gives the analytical continuation of the t-channel unitarity at large values 
of energy. For the BFKL Pomeron exchange we have 

Njp{Ko,ko,b;Y;Yo) = Njp{L,b,Y) = j ^ njp (7, 5) e'^^'^^^) ^ " (i"^) ^ (2.13) 

dH' [dl [ ^ [ ilnjp(7^6-6')e^«^(^)(^-^')-(l-^^)(^-')n^(7,60e^^^(^')^'-(l-^')' 



where L = In (^K^/kQ^ and / = In {k"^ /k'o). Kq and kQ are the momenta of the dipoles in the projectile 
and the target, respectively. Integrating over / we obtain that 7^ = 7. Considering y' = Yq we obtain that 

NHL,l,b;Y;Yo) = 27r^as(^HL - l,b,Y - Yq)) = N {b; L - l,Y - Yq) (2.14) 



Assuming that Eq. ( 2.13| ) and Eq. ( p.l4| ) hold in the general case but not only for the BFKL Pomeron 



exchange, we reduce Eq. ( 2.1C| ) to the following equation for the amplitudes 



+ 5s|2iV(L-i,(,,y'-r') (l + i)'^ "('.6,1") + (| + i)'^A'={i,(,,y') 1 

The second equation that stems from variation with respect to ^ {I, b, Y') has the same form as Eq. ( p.l5|) . 
In Ref. [5] we solve Eq. (|1|) in semi-classical approximation assuming that 



N{l,b,Y') = e^^''''"^') = e"('''''^')^'-(i-^('''''^'))' (2.16) 



and using the method of characteristics. In Eq. ( 2.16| ) we consider that u}{l,b,Y') = dS {l,b,Y') jdY' 



and 7 (^, 6, F') = dS {I, b, Y') /dl are smooth functions of Y' and I (see Ref. [5] for more details). We found 
that for any value of zy there exists the solution at large / which is very close to the solution of the linear 
BFKL equation. In particular this solution has a critical characteristic for 7 = 'jcr that can be found from 
the following equation [7, 18] 
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xilcr) _ xilcr) T^-j^Qj-Q xil) = 2^(1) — (7) ~ ^ (1 ~ 7) ^ kernel of the BFKL equation 



(2.17) 



1 - d-Jcr 

with ip{z) = dliiT(z)/dz and T(z) is the Euler gamma-function. 
The equation for the saturation scale looks as follows 

. = ln(^^^) = as^r - I (2.18) 

In the vicinity of the saturation scale but for 2; < the scattering amplitude shows the geometric scaling 
behaviour [19] i.e. it depends only on one variable (z) instead of three: Y', I and b. For the Balitsky- 
Kovchegov equation the geometric scaling behaviour of the scattering amplitude is the typical feature inside 
the saturation region (see Ref. [20,21]). In this paper we are going to solve Eq. ( 2.15| ) treating (z; Y') as 



a function of two variable: z (see Eq. ( |2.18| )) and Y'. The choice of the variable shows that we believe that 
the scattering amplitude inside the saturation region has the geometric scaling behaviour and the initial 
condition for this solution can be found from the solution outside of the saturation scale, namely, 

N{z = 0,Y') = No; ^^U=o = (1 " 7cr) iVo = 1 - e"^" (2.19) 

Recall that on the critical trajectory, the amplitude is constant and in the vicinity of the saturation scale it 
is proportional to A^o exp ^ (1 ~ 7cr) ■ However, introducing a dependence on Y' we are going to check 
whether the assumption on the scaling behavior of the amplitude is correct and within what accuracy. The 
initial conditions at y = we set using the McLerran - Venugopalan formula [9] (see term {. . . } below), 
namely 

N{z,Y' = 0) = - e-«^(^-0)-?2/4| = l-r{0,e/Ql{Y' = 0)) (2.20) 

2.2 Solution inside the saturation domain: general approach 

For finding the solution inside the saturation region we will use a method proposed in Ref. [21] (see also 
Refs. [5,6]): we introduce function (j){z) as follows 

N{z,Y') = \[ ^^'(1 - e-'^^"''^'^) + No (2.21) 

and assuming that function d(j) {z, Y') /dz is a smooth function we will find the solution to Eq. ( |2.15| ). 
The smoothness of function d(f){z,Y') /dz means that 



Using Eq. ( 2.22| ) as well as the properties of the BFKL equation (see Ref. [5] for details) we obtain the 



following equation 
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w + £ (7) - (7) e"'^^^^ = 2N{z,Y') + 2N^z,Y') = 2N {z,Y') + 2N^ZY - z,Y - ¥') (2.23) 



Where 



7 = ^ and w = _ ; 2.24 

oz as oY' 

(l){z = 0,Y') = {l-jcr)No = 0.63 iVo and (l){z,Y' = 0) = ^oe^; (2.25) 

\ Xilcr) 1 1 + 37 

1 - 7cr 7 7(1 + 7 



1 + 67 + 77^ ^ 

•^(^) = 7(1+7)2 ' (2.27) 



Introducing N = 2N we can rewrite Eq. ( p.23| ) in the form 
N'Jy, {z, Y') 



{l-N'Az^Y')) 



(2.28) 



c \ , ^-j"'"^^^ ^ ]-^\ f ""''-"'^'^ ^ \ i^-^'A^y)) - ^i^y) - N{zy-z,Y-r) 

[l-N'Jz,Y'))) \[l - N'Jz,Y')) J ^ ^ 
Initial conditions for this equation are given by Eq. (|2.25| ). 
2.3 Asymptotic solution 



Demonstrating the main features and the problems that we face solving Eq. ( 2.23| ) we, first, investigate 



the specific case: (j)(z,Y') increases at large z and it has a geometric scaling behaviour (p{z,Y') = (l){z). 
Using these assumptions we can simplify the general equation (see Eq. ( ^.23D ): 

£(7) = 2N{z) + 2N{zY-z) = ZY (2.29) 

Function £ (7 = d(l)/dz) is shown in Fig. Q (black line). One can see that we have the following asymptotic 
solutions [5]: 

1. For zy > z™™ ~ 6.3 there is a solution at (j)'^ {z) = 7 (z) <C 1(7^0); 

2. For all zy we have solutions (j)' {z) = nz n), where n = 1, 2, 3, ... ; 

3. For Zy > z^^^ ~ 6.3 we have the following solution: 

^{z) = ^-^P^ZYZ (2.30) 

X (7cr) 
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X(r) 




Figure 1: Function £ (7 = di/i (z) /dz) versus 7 (black line). Blue curve shows £ (7) + tt cot (7r7) — I/7 while the red 
one shows £ (7) — zy ■ 

In Fig. m the blue solid line shows 

t = £(7) +7rcot(7r7) - 1/7. (2.31) 

All poles at 7 = n are excluded in this function while the behaviour at large 7 this function has the same 
as C (7). This asymptotic behaviour at 7 close to zy/C (7cr) can be translated in Eq. ( |2.30| ) (see also below 



Eq. (pD and Eq. Q). 

Therefore, one can see that we face two major problems in searching the solution: (i) at first sight we 
have infinite number of solutions even in this simplified case; and (ii) we need to find the solution if we 
exclude the singularities in £ (7) . 

The infinite number of solutions contradicts the common sense intuition that the physical problem has 
the only one solution being formulated correctly. In the next section 3 we show mathematical arguments 
that will discriminate different solutions and will select the only one solution. This physical solution will 
be found in section 4. 

3. Semiclassical solution for (p (z, Y') 
3.1 Equations 

For large cp {z, Y') Eq. ( p. 23] ) can be re-written in the form* 

lD + £(7) - .F(7)e-* - ZY = + g {-/,(!))- ZY = (3.1) 



*We will denote below by zy the sum zy + 4 7Vo and hope it will not cause any difficulties in understanding. 
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Figure 2: Two kinematic regions: ili and 1^2 • 
The solid lines describe two sets of trajectories 
in fli and fl2- The trajectories in Qi start at 
z — with given value of 4>o- The trajectories 
in start at any point of z at y = 0. The 
initial value of 70 — (t>o exp(zo)- The dotted 
red line is the common characteristic line for 
both fii and ^l2- It has the form (see Eq. (4.5) 
and Eq. (0) below): Y' = z/Ci-jcr) ■ 



We solve this equation in semi-classical approximation as- 
suming that (z, Y') is a smooth function of both variables z 
and Y'. It is known (see Refs. [22] and references therein) that 
for the equation in the form 

F{Y',z,c^,-f,Cj) = (3.2) 

with smooth (j). We can introduce the set of characteristic lines 
: z{t), Y'{t), (j){t), u}{t), and 7(t) which are the functions of the 
variable t ( artificial time), that satisfy the following equations: 

(i£ ^ ^ ^ dg (7, (p) 
dt ^ 
dY' 



(2.) 

(3.) 

(4.) 
(5.) 



dt 
d4> 
lit 

d'j 

'dt 

du; 

~dt 



Fa, 



1 



dg 



+ Cj 



u 

-{F, + jF^) = -7-^(7) e-^ 
-- -{Fy'+uF^) = -cD^ (7) 6-^(3.3) 
Y'. Taking the ration of Eq. (|3.c 



From Eq. ( |3.3D -2 one can see that we can introduce t 
Eq. ( |3.3D -5 we obtain that 

u{Y') = IC{jo)l{Y') 

where 70 = 7 {Y' = 0) is the initial value of 7 which has to be found from the initial conditions. 



-4 to 



(3.4) 



Plugging Eq. ( |3.4D into Eq. (3J) we reduce the system of equation (see Eq. (3^)) to the following set 
of the equations: 

(3.5) 
(3.6) 
(3.7) 



£(7) - ZY + 7 /C (70))/ J" (7); 
dz{Y')/dY' = £^(7) - 7-^(7) e-^ 
dj/dY' = -j{c{j) + lC (70) 7 - zy) = T (7, 7o) ; 



It turns out the general set of trajectories can be divided in two groups: the trajectories that start at 
z = and at arbitrary Y' ( the vertical axis in Fig. ^) and the trajectories which starting points lie on the 
horizontal axis in Fig. |2| at Y' and arbitrary z. We need to consider these two sets separately. 

3.2 Solutions in fii 

In i7i we need to use the initial condition </) (z = 0) = 0o to obtain the equation for K (70)- 

/C(7o) = -(- £(70) + -F(7o) e--^" + zy) (3.8) 
70 V / 
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Fig. l-a Fig. |-b Fig. l-c 

Figure 5: exp (—0 (7, 70)) versus 7 for 70 = . . . l(Fig. |5]-a), for 7 = 1 . . . 2( Fig. |]-b) and 7 = 2... 3(Fig. ||-c) in 
the region fli. iVo = 0.1 and zy = 200. 



Substituting Eq. (|3^ 
the region l^i. 



into Eq. (|3.5D -Eq. ( p. 7] ) we have the final set of equations for the trajectories in 



Tiyo, To) 



In Fig. ^ we see that function T(7o,7o) is neg- 
ative for all values of 70. It means that j {Y') falls 
down and becomes frozen at the value 70 = 70,1 (70) 
at which ^(70,70) = (see Fig. |^. One can see 
that for < 7 < 1 actually equation r(7o,7o) = 
has two solutions 70,1 and 70,2 but both are smaller 
than 7o. For 7 > 1 we have the only one solu- 
tion 7o < 7o. Actually, 70 is very close to 70. In 
other words, 7 decreases very fast to 70. At 7 = 70 
T (7 = 7o, 7o) = 0. Since the scattering amplitude 
should be less that unity from the s-channel unitar- 
ity we expect that exp (—0 (z, Y')) in entire kinematic 
region should be positive and not equal to zero. Fig. ^ shows that we have a different behaviour of the 
solution for exp (^—(j)(^z {Y')^ ) on the trajectories, which contradicts our expectation from the physics point 
of view. 




Figure 3: r(7o,7o) versus 70 in the region ili 
0.1 and ZY = 200. 



An = 



Therefore, we can conclude that the semi-classical solution that satisfies the physical criterion, does 
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T(y) 




T(y, 1.5) 
T(y. 1.8) 
T{y, 1.98) 




T(y, 2.5) 
7(7, 2.8) 
7(7, 2.98) 



Fig. 



Fig. 0-b 



Fig. §c 



Figure 7: T (7, 70) versus 7 for 70 = . . . l(Fig. ^a), for 7 = 1 . . . 2( Fig. |^b) and 7 = 2... 3(Fig. ^c) in the region 
^2. TVo = 0.1 and zy = 200. 



not exist at any 70 in the region l^i when 7 close to n with n = 0, 1, 2, 3, . . . . 
3.3 Solutions in CI2 

In the region Q2 we have to repeat our analysis since the equation for /C (70) should be based on the 
initial condition: </> = zq, Y' = 0) = <j)Q exp (^^o)- Since 70 = (p{z = zq, Y' = 0) for this initial condition 
Eq. (|3.8D can be re-written in the form 



^ (70) 



-(- £(70) + J- (70) e-^o + zy 
70 ^ 



(3.9) 



Tiro, To) 



70 



One can see in Fig. |6| that T(7o,7o) in this re- 
gion is negative. Therefore, in this region as in region 
Oi 7 (z) falls down on the trajectory for all trajecto- 
ries. It turns out that for each value of 70 function 
^(7) 70) vanishes at 70 = 7o,i (70) (7^(70,70) = 0) 
(see Fig. |7|). Therefore, the solution in region Q2 has 
the same properties as the solution in region ili (see 
Fig. 1^ and Fig. |. 

Thus we can repeat the same conclusions for the 
region 0,2 as for Qi: there is no solutions that satisfy 
the physical criteria. 

4. Asymptotic solution with 7 y Zy j {^s j 



Figure 6: T(7o,7o) versus 70 in the region Q.^- Nq 
0.1 and zy = 200. 



Finally, the only solution which we need to consider is the solution with large 7 = d(j){z,Y') /dz 
zy/C (7cr) ^ 1- Eq. ( ^.23 ) reduces to the simple form 



^iz,Y') 



+ C(7cr 



dasY' dz 
with the following initial conditions: 

^{z = 0,y') = (/-o; 



zy 



where 



X ilcr) 

1 - 7c 



{z,Y = 0') 



(4.1) 



(4.2) 
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70=0.5 
70=0.8 
70=0.98 




70=1.5 
70=1.8 
70=1.98 




70=2.5 
y0=2.8 
70=2.98 



Fig. §c 



Fig. 



Fig. §c 



Figure 8: exp (—0 (7, 70)) versus 7 for 70 = . . . l(Fig. |8|-a), for 7 = 1 . . . 2( Fig. ||-b) and 7 = 2... 3(Fig. ||-c) in 
the region il2. iVo = 0.1 and zy — 200. 

We need to consider two separate kinematic regions Vl\ and for searching the solutions for this 
equation (see Fig. |2|): Y' > z/C (7cr) and Y' < z/C (7cr) [23]. In both regions the general solution to 
Eq. takes the form 

^{z,Y') =\zY {Y' + z/C{-fcr)) + G{Y' - z/C{-fcr)) (4.3) 
where G is the arbitrary function that has to be found from Eq. (|4.2| ). In region Q.2 one can see that 

G2 = (/^o exp ( - C (7e.) [Y' - z/C (7^.)) ) (4.4) 



from the second of Eq. ( [4. 21 ). However, in region r^i we need to use the first of Eq. ( |4.2| ) and we obtain 
that 



1 



ZY {Y' - Z/C(7cr.)) + 



(4.5) 



One can see that at Y' = z/C {jcr) two 0's: (pi (z) = zy z + 4>q and <j)2 (z, Y') = ^ zy {Y' + z/C (7cr)) + 
(po exp G (7cr) (X' — z/C (7cr)) ^ are equal, providing the needed matching. 

Therefore, the simplified Eq. ( [4.1D leads to the geometric scaling solution for Y' > z/C {'jcr) while 
for Y' < z/C{-f cr) we have a solution with explicit scaling violation. 

Armed with the asymptotic solution given by Eq. ( [4.5D and Eq. (^^) we study the numerical solution 
to Eq. (2.23) in region Qi assuming the geometric scaling behaviour and considering the following iterative 
procedure. Plugging in Eq. ( 2.23| ) (j){z,Y') = 4>{z) and Eq. (2.31) we see that this equation takes the 
following form in this region 



X(7c 



{z) 



1 - 7c 



dz 



2 7s + In 



dz 



zy 



H{z)-H{zy-z) with H{z) 



dz e 



/^-</,(i-i)(y) 



(4.6) 

where 7^; is the Euler constant (7^; = 0.5777216) and where c/)^*^ is the solution of the i-th iteration of the 
equation. Deriving Eq. ( [4. 61 ) we used that C 2^^^ 7^ + 111(7). 

We solve Eq. ( |4.6D using iteration procedure with 

',(*=o) = / (^^ + 2 (7s + Inzy)) (1 - 7,,)/x {-far) } z. (4.7) 



- 11 - 




N'.fe) 



T(z) 




Fig. U-h 



Fig. §c 



Figure 9: The third iteration of Eq. {L6): solutions for functions (f) (z)(Fig. |g|-a) and the scattering ampUtude N'^ (z) 
(Fig. ||-b) versus z. In Fig. ^c we plot function T (z) (see Eq. (4^)). In the figure the following parameters are 
taken: = 0.1 and zy = 200. 

The result of the third iteration is shown in Fig. One can see from Fig. ^-a that the solution is very 
close to (f)^^^ given by Eq. (O). N'^ reaches the unitarity limit at small values of z{z < 0.5). In Fig. ^c we 
plotted the difference between the l.h.s. and the r.h.s. of Eq. ( 2.23| ). This difference can be written in the 
form: 



... a1 l\f{i) ( y\ 

where e""^" = 1 - diV» (z) /dz and d^^'^/dz = ^ ' 



^^2 -/{l-dM^{z)/dz) 

T^*) (z) turns out to be small ( less than 1) leading to the accuracy of the solution about 2%. 

For the region 1^2 (see Fig. ^) we solve a more general equation for (/>(z,y). In this case Eq. ( |4.6D 



takes the form 

d^{z,Y') , Xilcr) d^{z,Y') 



dasY' 



+ 



1 - 7c 



dz 



2 7ij + In 



dz 



= ZY +4No- H{z,Y')-H{zY-z,Y-Y') 

(4.9) 



where H{z,Y') = dz'e '?^(^''^'). The initial conditions are given by Eq. ([4.2|). 

The numerical solution is shown in Fig. One can see that at large values of Y' solution approaches 
the geometric scaling solution of Eq. ( ^4.61) . We can use this physical solution, shown in Fig. [l^ and Fig. 
to study the N-N scattering in the next section. 

5. Nucleus-nucleus scattering amplitude 

We need to specify term Se in Eq. (2J) for calculating the scattering amplitude with a nucleus. This term 
determines the interaction of the BFKL Pomeron with the nucleons of the nucleus and it has been written 
in Ref. [1] in the following form 

Se = j dHd'^k(^<^{k,b,Y' = 0)TAAk,b) + {k,b,Y' = Y) TAAk,b)^ (5.1) 
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2000 




Figure 10: Function <l){z^Y') at different values of Y' as function of z. In the figure the following parameters are 
taken: iVo = 0.1 and zy = 200 



where 

/+00 p 
dzp{z,b) and a (k) = / n{k,h)d% (5.2) 
-oo J 

In Eq. ( |5.2| ) p is the density of the nucleons in the nucleus and a (k) is the cross section (imaginary part 
of the forward scattering amplitude) of the dipole with the nucleon at low energy while n{k,b) is the 
dipole-nucleon scattering amplitude. 

As we have eluded, in our treatment of nucleus-nucleus interaction we use the Glauber-type approach 
[26] integrating over all impact parameters of dipole-dipole and dipole-nucleon interaction since they are 
assumed to be much smaller than the impact parameters of nucleon-nucleon scattering (see for example 
Ref. [27] in which this approach is discussed in details). The latter is of the order of Ra_, which is larger 
than the nucleon radius and the sizes of all interacting dipoles. As we have discussed, we assume that 
instead of the real nuclei we are dealing with the nuclei that consist of mesons made of heavy quark and 
antiquarks. For such mesons we can calculate n {k) in the Born Approximation of perturbative QCD. In 
coordinate space n (r) = (27ra|Ci?/A'c) r^ln (i?^/r^) Q {R — r) where R is the radius of the meson which 
is of the order of 1/ {as (Mq) Mq) where Mq is the mass of heavy quark. In the momentum space we have 
the following ta {k, b) 

rMKb) = m . (*) = m ^ ^^^ip^ (5.3) 

r determine the initial conditions for the amplitude (A;, b, Y') and A^ (k, b, Y'), namely 



N {k, b, Y' = 0) = TA, (6, k) ; At {k, b, Y' = Y) = ta, {b, k) ; 



(5.4) 




(-1/3!) 0/2Q/2 Q/2 

Figure 11: The example of two nuclei reducible diagram in the scattering amplitude for nucleus- nucleus interaction. 

For simplicity, we consider the cylindric nuclei for which the h dependence is given by {Ra — b). In 
this model 

SAib) = 2poRA&iRA-b) (5.5) 

In this simple model the entire b dependence of (fe, b, Y') and A^^ {k, b, Y') turns out to be the same 
as in initial condition leading to 

N{k,b,Y') = SA{b)a{k,Y') = SA{b) {k,Y') 



N^k,b.Y') = SA{b)a^k,Y') = SA{b) -^M^k,Y') 



(5.6) 



where kg is the typical transverse momentum in the proton (fco oc asMg). Functions Af and A/"^ are 
dimensionless and for them we have the geometric scaling behaviour inside the saturation region [20] and 
in the vicinity of the saturation outside of the saturation region [19]. 

Using p {z = 0,b = 0) = po = 0.171/ fm^ and Ra = 1.2 A^^^ fm we can estimate the value of r at 
= for the gold: r = 8po RA/i^n M^) = 0.32(1/ fm'^)/M^ . One can see that A^o = 0.1 which we used 
for the numerical solution can be reached even at A; = if Mq ~ 0.36 GeV. 

Calculating the scattering amplitude, we need first to sum all two nuclei reducible diagrams (see Fig. 
for example). In such diagrams we can single out one or more states with two nuclei in the s-channel (see, 
for example, two such states in Fig. 0). 

This sum can be written as follows 



ImA (b, Y) 



1 



e 2 



in{6,y) 



(5.7) 



The notation : Q in Eq. (^.7|), is introduced to be opacity for the nucleus- nucleus scattering in the Glauber- 
Gribov approach [26]. 



Using Eq. ( ^.l]) one can see that the equation for the opacity ft has the following form (see Fig. 12) 
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Figure 12: The equation for opacity fl. 



n 



(6, y) = J dH' cPk TA^ (b - b', N {h' , kX = Y) = j dH' cPk (b - b', k, Y' = o) ta^ {b\ k) 

(5.8) 

Plugging in Eq. (p.6p we can rewrite Eq. (|5.8p in the form 



VI (6, Y) = Taa (b) J (fk a(b - b' , k^ a {k, Y' = Y) = Taa {b) ^ j (fka{b - b' , k^ N [k, Y' = Y) 

(5.9) 



where 



Taa (b) = {2poRAf j dH'Q {b' - Ra) 9 (|6 - b'\ - Ra) (5.10) 



{2poRAf Ra I 2 arccos ( — - ) - ;r^Wl 



b \ b b^ 



2RaJ 2RaV AR 



2 



As we have discussed the solution for M [z) at z < Zmin and M'^ {Zy — z) for zy — z < Zmin are the 
solution to the linear BFKL equation which can be written in the following form 

0<z<z^in N{k,Y) = ^N{z) = J (fk'a{k') aBFKL{k/k',Y' = Y) (5.11) 

alCp 1 f dk'^ 



< ZY-z < Zmin {k^ = {z) = j (fk'a{k') 4^^^ {k/k',Y' = Y) (5.12) 



In the vicinity of the saturation scale Mbfkl and J^bfkl -^1- ( 5.11 ) and Eq. ( 5.12 ) can be written 
in the simple form [20]: 



Mbfkl {k / k', Y' = Y) = Noe^^~''^^^'' and Af^pj^^ {k/k' ,Y' = Y) = N^e^^-^^^^i'^-^'^ (5.13) 
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N 



Y' 




Figure 13: Graphic form of Eq. (5.16). Solid black lines denote the Pomerons and their interactions that contribute 
to N, while the same contributions for N'^ are shown in red. 



where 



1-7. 



(5.14) 



cr 



Taking the integration over k' in Eq. (5.11) and Eq. ( |5.12| ) we obtain: 



^ ^ ^ Zmin 
< ZY - Z < Zmin 



M{z) 
M^{z) 



a^Cp No 



^exp ( (1 - 7cr) z) 

1 -Jcr ^ ' 



a\CF No 

TtNc I- 'Jc 



exp ( (1 - 7cr.) {zY - z) 



(5.15) 



It should be noticed that in the above equations / is defined as / = In {ji^ / {asMq)"^^ (see Eq. ( |2.1^ 
and Eq. (|]l|)) with /eg = {asMqf . 

Having these equation in mind we use a generaUzation of Eq. (^.9D 



(5.16) 



a (b; Y) = Taa (b) j dl A/'t {Y -Y'-L-I) M {Y'- l) 



Eq. (5.16|) has been proposed and discussed in Ref. [17] and for the nucleus-nucleus scattering it is 



illustrated by Fig. 13. From this figure one can see that arbitrary BFKL Pomeron diagram for the case of 



nucleus- nucleus scattering can be written as the product of N^ N. The extra Pomeron contribution that 
could connect two sets of diagrams for N^ and A^(shown in black for N and in red for A^^ in Fig. [l^ ) 
leads to a small corrections (see Ref. [4] for details). Y' is the arbitrary rapidity which is chosen from the 
condition z = Zmin- 

Using Eq. ( ^^ for A/'"'' {Y — Y' ; L — I) we can reduce Eq. ( p.l6| ) to the form 



n{b,ZY) = TAA{h) 



alCp No 



vr 



N^ {asMqY 1 - lev Jzy+z^.n 



rf^e(i-7cO(^y-^)_V'(z) 



(5.17) 



This equation stems from the geometric scaling behaviour of the solution to the equations for A^ and 



In the kinematic region z < 2zmin both A^ and N"^ in Eq. (|5.16 ) are the solution of the BFKL equation 
and in this region can be written as follows using Eq. ( 2.12| ): 
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Zy < 2 



Taa jb) / a|CF y 
{asMgf \ttNJ 



(1 - Icr)' 



■exp((l -7cr)zy) 



(5.18) 



For large z > 2zjnin ^{z) is given by solution of Eq. ([4.(]|). For estimates we can use for N (z) in this 
region the solution in the form 



dz' 1 



-</.('=0) (z') 



+ J^(z„ 



(5.19) 



where (/)^* is given by Eq. ( [4.71) . The expression for reads as follows 



Zy > 2 z„ 



fi(6,Zy) 

iVo 



TTNr 



+ 



(1 - Icr) 



{asMgY 

exp (^2 (1 - 7cr) Zmin^ | 



(iz'f 1 



-(/)(^=0)(2') 



(5.20) 



The result of our estimates using Eq. ( p. 18 ) - Eq. ( p. 201 ) is plotted in Fig. |T^. One can see that starting 
with small values of zy > 8 the amplitude at 6 = is close to the unitarity bound (see red line in Fig. [Tl|) . 
For comparison we plot in Fig. 14 also the amplitude which corresponds to the the exchange by one BFKL 
Pomeron in nucleon-nucleon scattering: 



Ap {AA;zY,b) 



1 — exp ( — -fip (zy , 6) 



(5.21) 



with Op (zy, b) given by Eq. ( |5.11|) at any value of 



zy. 



The main conclusions that we can make from Fig. |lj is the fact that the exact solution for nucleus- 
nucleus scattering shows the same corrections as Glauber- Gribov formula. This happens because at zy = 
2zmm the amplitude turns out to be very close to the unitarity bound. However, it should be stress that 
actually the exact solution leads to slower approaching to the unitarity bound since it turns out (see Fig. [Tsl ) 
that Q (zy , b) of Eq. ( |5.17| ) is much smaller than Qp (zy, 6) of Eq. ( 5.1^ ) that corresponds to the exchange 
of one BFKL Pomeron at the entire kinematic region. 



6. Conclusions 



The main result of this paper is the solution in the form of Eq. ( p. 21 ) for amplitude (Y, Y', I) = N (z, Y' 



with (j) (z, Y') given by Eq. ( [4. 3D and Eq. ( |4.6| ) and with the initial function (p^ ' (z, Y') determined by 
Eq. ( [4.2[ ). We obtained this solution in two steps. First, we assume that Eq. ( p. 14]) {N {l,b;Y;YQ) = 
{b; L — l,Y — Yq)) which is correct for the linear equation, is valid for the solution of the non-linear one. 
This conjecture allows us to reduce the system of two equations to one functional differential equation. 
Second, we solve this equation. The semi-classical approach was used to select out the infinite number of 
parasite solutions and the final solution was found analytically at large values of Y' as well as numerically 
in the entire region of Y'. 
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Im A(zy, b) 




Figure 14: The scattering amplitude for gold-gold interaction versus z = zy at different values of b: the solid black 



line line is the result of this paper (see Eq. (5.20)), the red dotted curve describes the Glauber-Gribov formula for 
one BFKL Pomeron exchange (see Eq. (5.21)); and the red line shows the unitarity bound. The curves correspond 
6 = 0,5 fm and 10 fm going from left to right, 



n(zy, 0) 



n(zy, 5) 





10) 




Figure 15: The opacity f](zy,6) of Eq. ( 5.2C ) and ftp (zy) of Eq. ( 5.11 ) for gold-gold interaction at different b 
versus z = zy- black line is the result of this paper (see Eq. (5.18)), the red dotted curve describes the contribution 
to D, from exchange of one BFKL Pomeron. The values of b are given in fermi. 



The solution, that has been found, looks as being different from the numerical solutions found in 
Refs. [2,3], but it shares with them several common features. In particular, this solution exists only and 
z < Zjnin ~ 4. In our solution the amplitude A^^ {Y,Y',l) = N (zy — z). We have no proof that this 
solution is unique but we believe that the physics problem has the only one solution. 

Using this solution we found the nucleus-nucleus scattering amplitude as a function of energy for 
'theoretical' nuclei. This 'theoretical' nucleus consists of the dipoles that are made of heavy quarks and 
antiquarks. We can use the perturbative QCD treating this nucleus. In reconstruction of the scattering 
amplitude we use Eq. ( [5.16| ) (see Fig. |l2|) which is the form of t-channel unitarity constraints and follows 
from Eq. ( ^.1| ) for Se term in the action of Eq. (p.l|). 

We showed that for zy < '^Zmin the exact calculation give the same result as the Glauber-Gribov 
formula which is widely used for nucleus-nucleus scattering. It turns out that z^in ~ 4 is so large that our 
estimates lead to the result close to the Glauber-Gribov formula in the entire kinematic region of energies. 
Only at large impact parameters the exact formula for nucleus-nucleus scattering gives visible deviations 
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from the estimates based on Glauber-Gribov approach. We beheve that this observation is important for 
all practical estimates based on Glauber-Gribov approach. 

In spite of the fact that the solved problem is still far away from the real physics environment we hope 
that our solution gives a reasonable first approximation to approach the nucleus-nucleus scattering for the 
nuclei that exist in reality. 
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